H band Image of a Planetary Companion around HR 8799 in 

2002 1 



Misato Fukagawa^, Yoichi Itoh^, Motohide Tamura'^'^, Yumiko Oasa^, Saeko S. Hayashi^'^, 
Yutaka Fujita^, Hiroshi Shibai^ and Masahiko Hayashi^'^ 

ABSTRACT 



The discovery of three planetary companions around HR 8799 (IMarois et al. 



20081 ) marked a significant epoch in direct imaging of extrasolar planets. Given 
the importance of this system, we re-analyzed H band images of HR 8799 ob- 
tained with the Subaru 36-elements adaptive optics (AO) in July 2002. The low- 
order AO imaging combined with the classical PSF-subtraction methods even 
revealed the extrasolar planet, HR 8799b. Our observations in 2002 confirmed 
that it has been orbiting HR 8799 in a counter-clockwise direction. The fiux of 
HR 8799b was consistent with those in the later epochs within th e uncertainty 



of 0.2 5 mag, further supporting the planetary mass estimate by iMarois et al. 



(]2008f ). 



Subject headings: stars: individual (HR 8799) — planetary systems 



1. Introduction 

The number of extrasolar planets discovered is increasing but primarily with indirect 
methods such as measurements of radial velocity attributed to the wobble of parent stars 
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have been achieved with Spitzer Space Telescope that cap tured the thermal infrared light 
during their secondary transits (jCharbonneau et al.ll2005l ). Direct imaging, in which the 
light from planets can be spatially separated from the host star, has long been desired 
since it opens up detailed characterization of the planets. The planetary-mass companion 
candidates have indeed been imaged aro und several young pre-main- or main-sequence stars 
( iLagrange et al.ll2009l : iKalas et al.ll2008l ). but the derivation of the planetary mass is often 
model-dependent or it is still uncertain if t he companion candida t e is gravitationally bo und 



to the host star (e.g., IChauvin et al.ll2004l : iNeuhauser et al.ll2005l : iLafreniere et al.ll2008l ). 



Most recently, iMarois et al.l (120081 ) reported one of the robust and exciting results: the 
discovery of a planetary system around an A5 star HR 8799 with direct imaging. The three 
Jupiter-like planets, HR 8799b, c, and d were detected at the separations relative to the 
star of 1'.'72, 0'.'96, and 0'.'62, respectively. Their multi-epoch astrometry from 2004 to 2008 
clearly shows that the companion candidates are not background objects but rather orbit- 
ing the central star in a counter-clockwise direction, probably moving in a nearly circular 
orbit. It is worth noting that their results in 2007 and 2008, that led to the significant de- 
tections of HR 8799c and d, were obtained using the angular differential imaging technique 
that effectively suppressed the sp ec kle noise . The host st ar , HR 8799 (HD 218396), has 
been recognized as a Vega-like star (jSylvester et al.lll996l ). Although the infrared exces s is 



van 



LeeuwenI 2007) is 



not quite large (Lir/L^, = 2.0 x 10~^), its proximity {d = 39.4 pc 
suited to resolve the structure in the inner part, and hence for direct imaging. Therefore, 
we coronagraphically observed it in 2002 with the Subaru Telescop e, originally intende d to 
find scattered light from the debris disk. Following the discovery by lMarois et al.l (120081 ). we 
revisited our Subaru data and applied thorough analysis suited for point source detection, 
since the astrometry in other epochs would have an importance for future accurate deter- 
minations of the planetary orbits, as well as for independent confirmation of the planetary 
companions. In the present Letter, we report the detection of HR 8799b in 2002. 



2. Observations 



HR 8799 {H = 5.28) was observed on 2002 July 19 (UT) with the near- infrared coron- 
agraphic camera CIAO (Coronagraphic Imager with Adaptive Optics; iTamura et al.ll2000l ) 
mounted on the Subaru Telescope. The image s were obtained in t he H band (A = 1.65 fim) 
with the Subaru 36-elements adaptive optics (ITakami et al.ll2004j ). The 0'.'5 diameter trans- 
missive coronagraphic mask and the 80% diameter Lyot stop without spider vanes were used. 
The sky was clear and the seeing was relatively good (~ 0'.'6 at optical) and stable except 
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for at the beginning of the observations. No special technique optimized for planet detec- 
tion was applied such as angular or spectral differential imaging. CIAO is equipped with a 
1024x1024 InSb ALADDIN II array. Since our observations were aimed to detect faint ex- 
tended emission from the debris disk, the medium-resolution mode (21 mas pixel"^) was used 
to obtain the better sensitivity. Nevertheless, the pixel scale of 21 mas was sufficiently small 
to measure the positional difference for the companion candidates aro und HR 8799 give n 
their large orbital motions (25 mas/yr for the most distant companion; iMarois et al.ll2008l ). 
The pixel scale and orientation of the detector that we adopt here are 21'.'2±0'.'4 mas pixel ~^ 
and 5°9 ib0?5, respective ly, derived from the observations in 2001 and 2002 of the Trapezium 
cluster (jitoh et al.ll2005l ). M15, and objects that have other nearby point sources. The dis- 
tortion is negligible and was not corrected. As a PSF reference star, HD 218172 {H = 5.93, 
separated by 1° from HD 218396 on the projected sky) was also imaged with the same occult- 
ing mask of 0'.'5 just before and after observing HR 8799. The AO guide stars were HR 8799 
and the reference star themselves, but they were too bright for the AO and their fluxes did 
not match in R band where the wavefront sensing was performed. Therefore, the neutral 
density filters were inserted in front of the lenslet array in the AO unit, and by selecting the 
appropriate number of the filters for each star, the similar performance in the atmospheric 
correction was achieved for both sources. The 5 exposures of 1 sec were coadded into one 
data frame, and 264 and 200 frames were obtained in total for HR 8799 and the PSF ref- 
erence star, respectively. The position angle of the spider pattern was varied with time in 
the images, since the instrument rotator was used with the altitude-azimuth telescope. HR 
8799 and the PSF reference star were not observed at the meridian (near zenith), and the 
spider pattern was rotated only by 4 degrees during the observation of HR 8799. The stars, 
which were tracked at optical, were gradually misaligned with the mask center with time, 
and the positional adjustment was made with the AO tip-tilt mirror. The focal plane mask 
was moved by 1" and the star was placed on the mask using the AO tip-tilt after taking 132 
frames for HR 8799 (dithering). FS 155 was observed on the same night as a fiux calibrator, 
without the occulting mask and without the AO loop closed to minimize the observational 
overhead. 



3. Data Reductions 

The obtained images were calibrated using IRAF for dark subtraction, fiat-fielding with 
dome-fiats, bad pixel substitution, and sky subtraction. The dark exposures were obtained 
before the observations were started, and the dome fiats were taken using an incandescent 
lamp at the end of the night. The constant sky level was estimated for each frame measuring 
the median value at large separations (> 5" from the star). 
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In order to extract the extended emission or companion candidates in the vicinity of 
the bright star, the PSF-subtraction was carried out. First, the careful selection of the 
usable data was required to reduce the effect of the temporal variations in the PSFs. As a 
measure of the PSF quality, we used the brightness ratio of the halo at r = I'.'S and the peak 
intensity, then eliminated the images with the ratio larger than 0.02 for both HR 8799 and 
the PSF reference. The 220 frames were selected for HR 8799, corresponding to the total 
integration time of 18.3 minutes, and the 141 frames (11.8 minutes) were adopted for the 
reference star. Second, the stellar position was estimated by measuring the centroid of the 
transmitted stellar light in the mask, and the stars were registered to a common position. 
Each four consecutive frames were averaged to improve the signal-to-noise ratio (SNR). 

We applied the following procedures in order to remove the stellar halo. 

1. Subtraction of the image rotated by 180°. 

The rotated image for HR 8799 was subtracted from the original unrotated one, but 
there remained an asymmetric positive and negative pattern in every frame since the 
PSF was not completely circular. In order not to cause the asymmetry, we had to 
slightly shift the rotated image. The amount of the shifts depended on each frame, 
but the maximum offsets are Aa = 0'.'065 and AS = 0'.'028. The relatively large shifts 
do not indicate that the centroid measurement of the stellar position was incorrect, 
since the spider pattern was clearly misaligned in the subtracted result with the re- 
registered frame. However, we conservatively included the averaged offsets, Aa = 0'.'021 
and A6 = O'.'OOOS, in the error of the astrometry for HR 8799b described in the later 
section. The subtraction was made for each frame of HR 8799, and the resultant 
frames were median combined (the panel (a) of Figure [T]). Furthermore, the same 
subtractions using the rotated images were performed for the PSF reference star. The 
result was scaled to the brightness of HR 8799, rotated to cancel the difference of the 
instrumental rotator angle, then subtracted from the result for HR 8799 to suppress 
the speckle noise (the panel (c) of Figured]). 

2. Subtraction of the PSF constructed from the reference star images. 

Five PSFs were built by averaging all the reference images, images obtained before or 
after observing HR 8799, images that have the brightness ratio of the halo (r = I'.'S) 
and the peak of < 0.01, and those with the ratio of < 0.008. The PSF subtraction 
was performed for each frame of HR 8799 using the averaged reference PSF with the 
higher SNR. The PSF was rotated to match the position angle of the spider pattern 
for HR 8799, and scaled to fit the halo brightness of HR 8799 at r = I'.'S. It was 
necessary to manually tweak the scaling factor and the position of the PSF to improve 
the subtraction. The best PSF was chosen so that it gave the smallest residuals. The 
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PSF-subtracted images were median combined and shown in the panel (a) of Figure [2j 
The real point source should be present in all the data frames. In order to extract 
such emission, the standard deviation of the pixel values was calculated at each pixel 
position (a noise map). The PSF-subtracted image was divided by the noise map to 
emphasize the components that were temporally stable (the panel (b) of Figure [2]). 

3. Subtraction of the median-filtered images. 

This procedure was applied in order to remove the emission with lower spatial frequency 
than point sources. The accurate measurement of the FWHM for a point source was 
impossible since there was no bright, unocculted star in the small field of view of 
CIAO. However, the light from the PSF core for HR 8799 was transmitted through the 
occulting mask, and the derived FWHM ((X'O? ± O'.'Ol) was a reasonable value for the 
A036 imaging using a bright guide star under good seeing condition. The consistency 
of the FWHMs between unmasked companions and a masked star was also confirmed 
for HD 200775 observed on the same night, but under worse seeing condition, i.e., 
the FWHM of O'.'l and the brightness ratio of the halo and the peak of >0.02. In 
order to extract the emission that has the similar spatial scale as a point source (0'.'07 
= 3.3 pixels), the median- filtered image was constructed using a median window of 
21x21 pixels and subtracted. The resultant images were median combined to reduce 
the peaky noise pixels (the panel (a) of Figure [3]). The reference frames used for the 
reduction 2 (the subtraction of the PSF of the reference star) were also processed in the 
same way, and subtracted from the filtered image for HR 8799. The each subtracted 
images were finally median combined and displayed in Figure [31 

4. Subtraction using the Locally Optimized Combination of Images (LOCI) algorithm. 



Lafreniere et al.l (120071 ) proposed the effective PSF subtraction algorithm in which ref- 



erence images are linearly combined with coefficients of the combinat ion estimated to 



minimize the residuals in each subsection of the image. Following iLafreniere et al. 



( 120071 ). the squared residuals of the subtraction of the reference PSF (O^) from the 
target (O"^) are expressed as = Y.ii.Oj - OfY = Yj^OJ - Y^k^^O^Y, where i 
denotes a pixel in the subsection and k is the index of the reference image. The coef- 
ficients for the fcth image, c*^, can be calculated from linear equations, Ac = b where 
^jk = Ei OjOf, bj = Yi OjOf. The PSFs for the reduction 2 (rotated and scaled 
ones) were used as reference images. Before applying the LOCI, we reduced the low- 
frequency emission that produced the huge mismatch of the PSFs, subtracting the 
median filtered image using a 21x21 box to calculate the median value. The sub- 
traction was made for each HR 8799 image in the region of interest (r > I'.'S), with 
the radial and azimuthal width of the subsection set to be O'.'S and 20°, respectively 
(the area corresponds to 79 PSF cores with 0'.'07 FWHM). The median combined re- 
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suit is shown in the panel (c) of Figure [2j We tried to apply the subtractions using 
smaller subsections for the inner region, but did not produce a significant detection of 
HR 8799c. 



4. Results 

The observations were originally designed to detect the light scattered from the disk. 
However, no extended emission was found beyond the boundary where the subtraction resid- 
ual started to dominate (r < 1'.'4). The detection limit is 16 mag arcsec"^ which was es- 
timated by placing an artificial disk having a uniform intensity of 2a with la noise in the 
PSF-subtracted image. 

Figure [H [21 and [3] show that ther e is a point-like emi ssion (emission B) very close to 



the location of HR 8799b reported by iMarois et al.l (120081 ). The emission is faint, but we 



confirmed the detection of the emission B based on the four different subtraction methods 
described in the previous section. The subtractions of the self-rotated images and self- 
median-filtered images also demonstrate that the detection is not an artifact by the mismatch 
of the PSFs of HR 8799 and the reference star (see the left panels of Figure [1] and [3]) . In 
addition, the emission B was found both before and after the dithering, indicating that the 
emission B is not a cluster of hot pixels (Figure [21 the bottom panels). The FWHM of the 
emission B was measured to be 0'.'07 ± 0'.'02, which is in good agreement with the FWHM 
estimated for the masked PSF core of HR 8799. Therefore, we conclude that the emission B 
is likely to be a real point source, and it corresponds to HR 8799b. 

The projected separation of HR 8799b from the parent star was measurable thanks 
to the transmissive coronagraphic mask. The centroid positions were measured, and their 
separation is obtained to be (Aa, A5) = (1'.'481 ± 0'.'023, 0'.'919 ± 0'.'017). The uncertainty 
includes the error of the plate scale and the uncertainty of the centroid measurements esti- 
mated based on the self-rotated image subtraction as described in the previous section. Since 
the centroid uncertainty derived in the DEC direction is unreasonably small (0.02 pixel), we 
further assumed the additional 0.5 pixel error. Those are the dominant error sources than 
the variation between the subtraction methods (PSF star subtraction, median-filtering, and 
the LOCI subtraction) and the frame selection (be fore or after the dith ering). The possibil- 



ity of a background object is clearly rejected as in iMarois et al.l (l2008l ): the proper motion 
of HR 8799 suggests that the relative separation of a background star to HR 8799 should 
be changed by 0'.'71 in 6 years, while it changed only 0'.'02 from July 2002 to 2008. The 
multi-epoch data are summarized in Table [H Our 2002 result also shows that HR 8799b was 
orbiting in counter-clockwise direction (Figure [1|) at a rate of 18t9^ "^^^ y^~^ from 2002 to 
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2004. The projected separation 1'.'743 ± 0'/029 corresponds to 68.7 AU assuming the distance 
of 39.4 pc. This is consistent with t hose in 2004-2 08 wi thin the uncertainty, indicating that 
the nearly circular orbit inferred by iMarois et al.l (120081 ) is plausible. 



The aperture photometry was made using the IRAF/APPHOT package with an aper- 
ture of 3 pixels in radius and the sky annulus of 5-10 pixels. The median-filtered images 
were not used for the photometry. The comparison with the photometry for HR 8799 may 
give a reasonable magnitude, but the transmission of the occulting mask has not been accu- 
rately determined, ranging from 1-2% based on the measurements during the commissioning 
observations, to 11% obtained with the laboratory experiment. Nevertheless, assuming the 
observed transmission of ~1%, the 3 pixel aperture photometry for the masked HR 8799 
suggests if = 18 mag for HR 8799b. Alternatively, the aperture correction from 3 to 6 pixels 
(approximately 4 HWHM) [AH = 0.26 mag) was applied, which was derived from the artifi- 
cial PSF modeled with a moffat function that fit the transmitted core of HR 8799. The flux 
uncertainty was estimated by considering the standard deviation of the photometric results 
for several data subsets as well as the uncertainty associated with the aperture correction, 
and the photometric errors for the flux calibrator and HR 8799b itself. In addition, the dif- 
ference of the photometric results using the reference star PSF and the self-rotated images 
was considered. T he inferred magnitud e is if = 18.07±0.25 m ag, which is in agreem ent with 
that measured by iMarois et al.l (120081 ). The agreement with iMarois et al.l (120081 ) suggests 
that the planet is n ot significantly variable in the interval of 6 years at its age (30-160 Myr; 
Marois et al.ll2008l ). In addition, the consistency further supports their derivation of the 
planetary mass. 

There cannot be seen other significant detections of point sources in the PSF-subtracted 
images due to the severe contamination by the subtraction residuals. We have to say that it 
is completely impossible to detect HR 8799d which is the closest one to the star (r = 0'.'6). 
There can be seen an emission at the expected location of HR 8799c (r ~ 1"), but it is 
strongly affected by the speckle noise that cannot be removed at such radii. HR 8799c is 
therefore not detected. 



5. Summary 



The imaging observations with the 36-elements AO rev ealed a planetary companion 
HR 8799b which has been discovered by IMarois et al.l (120081 ). Our Subaru data obtained 
in July 2002 gave the longest time coverage for the positional measurement when combined 
with the Keck results. It is confirmed that HR 8799b is bound to the host star and orbiting 
counter-clockwise. Note that these data points do not produce a smooth track (Figure S]) 
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and additional observations are required to derive the orbital solution. The photometric 
result is also consistent with the observations at later epochs. 
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Fig. 1. — HR 8799 in H band with the rotated image subtractions, (a) The median-combined 
image after subtracting the self-rotated images. The strong features running in north-south 
and east-west directions are caused due to the spider of the secondary mirror, and denoted 
with arrows, (b) The PSF reference star image with the subtraction of the rotated images, 
(c) The difference between the self-rotated subtractions for HR 8799 and the reference star. 
The field of view is 4" x 4". All the images are shown in the same display range. 
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Fig. 2. — HR 8799 in H band with the subtractions of the PSF constructed from the reference 
star images, (a) The median-combined PSF-subtracted image, (b) The PSF-subtracted 
image divided by the standard deviation of pixel values. The map was smoothed by the 
Gaussian function with its sigma of 1.6 pixel, (c) The median-combined images obtained 
before the dithering, and (d) after the dithering. The field of view is 4" x 4". The panel (a), 
(c), and (d) are displayed in the same range. The panel (b) was linearly shown in the range 
of 0.5-1.5, but if the Poisson noise component is subtracted, the SNR for the emission B is 
4. 
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Fig. 3. — (a) The subtractions of median-filtered self images, (b) The subtractions of the 
median-filtered reference images, (c) The result of the LOCI subtraction. The field of view 
is 4" X 4". All the images are shown in the same display range. 
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Fig. 4.— Or bital motion of HR 8799b from 2002 to 2008. The positions in 2004-2008 were 
measured by iMarois et al.l (120081 ) . 
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Table 1. Projected Separation for HR 8977b from the Parent Star 



date 




Aa n 


A5 n 


reference 


2002 Jul. 


19 


1.481 ± 0.023 


0.919 ± 0.017 


this work 




2004 Jul. 


14 


1.471 ± 0.005 


0.884 ± 0.005 


Marois et al. 


'2008) 


2007 Oct. 


25 


1.512 ± 0.005 


0.805 ± 0.005 


Marois et al. 


'2008) 


2008 Jul. 


11 


1.527 ± 0.004 


0.799 ± 0.004 


Marois et al. 


'20081 


2008 Aug 


12 


1.527 ± 0.002 


0.801 ± 0.002 


Marois et al. 


'20081 


2008 Sept 


. 18 


1.528 ± 0.003 


0.798 ± 0.003 


Marois et al. 


'20081 



Note. — The uncertainties of the plate scale and orientation 
for Keck/NIRC2 are well below the errors quoted in the Table 
jMarois et al.ll2008[ l. 



